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■ In the cold dark matter model of structure formation, galaxies are assembled hierarchically from 

' mergers and the accretion of subclumps. This process is expected to leave residual substructure in 

the Galactic dark halo, including partially disrupted clumps and their associated tidal debris. We 
develop a model for such halo substructure and study its implications for dark matter (WIMP and 
axion) detection experiments. We combine the Press-Schechter model for the distribution of halo 
subclump masses with N-body simulations of the evolution and disruption of individual clumps as 
they orbit through the evolving Galaxy to derive the probability that the Earth is passing through a 
subclump or stream of a given density. Our results suggest that it is likely that the local complement 
of dark matter particles includes a 1 — 5% contribution from a single clump. The implications for 
, dark matter detection experiments are significant, since the disrupted clump is composed of a 'cold' 

flow of high-velocity particles. We describe the distinctive features due to halo clumps that would 
be seen in the energy and angular spectra of detection experiments. The annual modulation of these 
features would have a different signature and phase from that for a smooth halo and, in principle, 
would allow one to discern the direction of motion of the clump relative to the Galactic center. 
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I. INTRODUCTION 



A cornerstone of modern cosmology is the hypothesis that most of the matter in the Universe is non-luminous and 
O I more than likely non-baryonic. The most convincing evidence for dark matter comes from spiral galaxies similar to 
our own. The rotation curve of a typical spiral is observed to be flat or rising slowly well beyond the Holmberg radius 
^ ' which encloses 90% of the galaxy's total luminosity. The usual interpretation is that spiral galaxies are embedded 
• • in extended halos of dark matter (e.g., The rotation curve probes the spherically-averaged integrated mass 

^ ' distribution - a flat rotation curve, for example, implies a density profile p{r) oc r~^. However, observations have so 
^> , far yielded relatively little information about the detailed structure of dark halos, e.g., their shapes, and even less is 
5^ ' known about their structure on small scales. 

pi , An understanding of the structure of the Galaxy's dark halo is paramount for dark matter detection experiments 
since they rely on halo model predictions to develop search strategies. Moreover, the analysis and interpretation of 
the experimental data, e.g., the derivation of particle physics parameter constraints from detections or non-detections, 
must contend with the large uncertainties in these models. The microlensing experiments designed to detect massive 
compact halo objects (MACHOs) provide a cautionary tale in this regard. It has been demonstrated that flattening, 
velocity space anisotropy, and rotation of the halo, all of which are consistent with other constraints, impact signif- 
icantly on the interpretation of the data ||^ . In addition, small-scale structure in the form of a clump or stream of 
MACHOs moving between us and the LMC could lead to an overestimate of the MACHO fraction in the Galaxy 

SilJ- 

One difficulty with the analysis of microlensing experiments is that there is no consensus on how a MACHO halo 
would form. The theoretical situation is significantly better for cold dark matter (CDM) candidates such as weakly 
interacting massive particles (WIMPs) and axions. In CDM models, theoretical arguments and N-body simulations 
indicate that structure forms hierarchically: small-scale objects collapse first and subsequently merge to form systems 
of increasing size. It is often assumed that once subsystems are assembled into galaxy-size objects, structure on 
subgalactic scales is erased through processes such as tidal stripping, violent relaxation, and phase mixing [Q. The 
resulting particle distribution of galaxy halos would be relatively smooth in both position and velocity space. This 
viewpoint is now being challenged by high-resolution N-body simulations of halo formation in cold dark matter models 
(e.g., 1^, 1^]): galaxy-size halos in these simulations have a great deal of surviving substructure in the form of subclumps 
and tidal streams. This appears to contrast with the observed structure of the Milky Way, which has relatively few 
luminous satellites. The discrepancy has stimulated the revival of models in which the dark matter has some property 
— e.g., interactions or non-negligible velocity dispersion — which prevents it from collapsing on small scales (e.g., 
pO| , pi] , p^). The alternative is that dark halos are in fact as clumpy as is seen in the simulations, but that gas 
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collapse and/or star formation have not occurred in a large fraction of the galactic subclumps and therefore they 
are not observed Moreover, there is mounting evidence for the existence of tidal streams at large galactocentric 
distances in the Milky Way ^4, |l5, as expected in a hierarchical cold dark matter model 



Theoretical predictions for WIMP and axion detection experiments have, for the most part, sidestepped the issue 
of halo substructure. The most widely used halo model for dark matter search experiments is an isothermal sphere 
with non-zero core radius and a truncated Maxwellian velocity distribution (e.g.,||l^). Although other halo models 
have been considered, in general they assume a smooth distribution function for the dark matter (e.g., jl^, 
An exception is the work of Sikivie et al. j2^, |2^, who assume that halo formation takes place through the 
accretion of successive spherical shells of matter. Spherical infall models have proved to be quite useful in the study 
of structure formation ||2^. For example, self-similar infall models ^ provide exact, fully non-linear solutions 
for the distribution function describing a collapsing system of collisionless particles in an expanding universe. A 
characteristic of these models is that the dark matter distribution at each location in the halo is composed of distinct 
streams of particles. The velocity dispersion within individual streams is small - a reflection of the fact that the initial 
distribution of particles is "cold" . Each stream is characterized by a different velocity that is directed either inward 
or outward, with ~ 60 individual streams at the position of the Sun. From a practical point of view, most of these 
can be viewed as components of a smooth background of particles whose velocity dispersion is given by the spread in 
velocities of the streams. However, streams associated with shells that have been accreted recently, i.e., those passing 
through the central regions of the Galaxy for the first, second, or third time, have speeds that are significantly higher 
than the characteristic 'thermal' speed of the bulk of the particles. They should therefore be observable as distinct 
features in the energy spectrum of dark matter particles detected on Earth. Some of the implications of this model 
for WIMP detection have been considered by Copi, et al. ||2^, and Gelmini & Gondolo | |2^ . 

The conclusions of Sikivie et al. |2j are valid only in so far as the spherical infall model actually applies to 

the formation of the Galactic halo. In this paper, we study the consequences for dark matter detection of substructure 
in the Galactic halo formed in the context of the more realistic hierarchical clustering scenario for cold dark matter. 
Halo substructure is assumed to be in the form of cold subclumps and their tidal debris (e.g., ^) - the remnants 
of hierarchical mergers and accretion - rather than cold infalling shells. While some of the observational signatures 
of substructure formed in the hierarchical model are qualitatively similar to those of the spherical infall model, the 
interpretation is quite different. At its core, hierarchical clustering is a random process. Our results are therefore 
expressed as a probability that a given feature will appear in the angular or energy distribution of events for a dark 
matter experiment. By contrast, the spherical infall model is deterministic once the profile of the initial density 
perturbation is specified. 

The most straightforward way to investigate this probability distribution would be to use standard N-body simu- 
lations and simply look at the phase space distribution in regions corresponding to the solar neighborhood. However, 
as illustrated by Moore [^, current simulations do not have sufficient resolution to fully probe the inner regions of 
dark halos. Therefore, an alternative approach is required. Our analysis combines results from phenomenological 
models of halo formation — namely the extended Press-Schechter formalism — with numerical simulations designed 
to follow the evolution of individual subclumps in a background gravitational potential. The ingredients of the model 
calculation are provided in Section II. Details about the gravitational potential assumed for the Galaxy and the 
initial structure of the clumps are given in Section III . Our numerical experiments are described in Section [V and the 
results are presented in Section 0. Experimental signatures for WIMP and axion detection experiments are discussed 
m Section 0. 

We note at the outset that our model for substructure in the Galaxy halo is simplified in a number of respects; 
nevertheless, we expect it to yield a qualitatively accurate picture of the possible deviations from smoothness in 
the local phase space of halo dark matter. For other recent investigations of the observational implications of halo 
substructure, see |3^, |3^, earlier work on possible implications of halo dumpiness for dark matter detection 
includes MM- 



II. PROBABILITY DISTRIBUTION OF HALO SUBSTRUCTURE 



We aim to estimate the probability that the local complement of dark matter particles includes a measurable 
contribution from a gravitationally bound clump or tidal stream. We focus on clumps that have made up to four 
orbits through the Galaxy by the present day — i.e., clumps that started to fall into the Milky Way at a redshift 
z < 1. The local distribution of dark matter particles is therefore divided into two components - a smooth background 
composed of particles that were accreted at early times (the substructure of which has since been erased by dynamical 
processes) and inhomogeneous material from recent accretion events. We assume that for z < 1 changes in the 
gravitational potential of the Galaxy are gradual and that clump-clump interactions can be ignored. Under these 
assumptions, recent accretion events can be studied numerically by evolving individual subclumps in a smooth, time- 
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dependent model potential. Our assumptions are based on a variety of arguments which indicate that the recent 
accretion rate onto the inner parts of the Galaxy has been relatively low. The coldness and thinness of the Galactic 
disk, for example, limit the infall rate of satellites since they can transfer energy to stars in the disk |3^. Measurements 
of the ages and metallicities of stars in the Milky Way's halo suggest that less than 10% of the halo stars come from 



recent merger events |38 . By contrast, numerical simulations and theoretical modeling imply that the mass of the 
extended dark halo has grown by a factor of 2 or more since z ~ 1. The conclusion is that most of the material 
accreted recently resides in the outer parts of the halo. Our analysis focuses on those few clumps that reach the inner 
regions of the Galaxy. 

Our results are expressed in terms of a probability distribution function dP/ dp, where dP is the probability that 
the density of dark matter particles in the solar neighborhood associated with a single clump or stream is between p 
and p + dp. In general, we can write 



dP 
dp 



7V(zta, / (^ta, rn,pj]p) dztadpjdm 



(2.1) 



where N' {ztg_, m,pj) dztudpjdm is the number of clumps accreted by the Galaxy with mass between m and m + dm, 
turnaround redshift between zta and zta + dztn and orbital parameters (e.g., angular momentum) between pj and 
Pj + dpj. The turnaround redshift is defined as the epoch at which a clump breaks away from the expansion and 
begins to fall in toward the center of the Galaxy. is modeled using the extended Press-Schechter formalism 
(see Section III). The dynamical evolution of accreted clumps within the Galactic halo is encoded in the function 
f{zta, rn,pj; p)dp, which gives the probability that a clump characterized by the parameters m, Zts., and pj contributes 
a density between p and p + dp to the present-day density of dark matter in the solar neighborhood. 

Before proceeding to the elements of the model, we first reduce the calculation of / to a more tractable problem. 
Consider a volume V representative of the solar neighborhood. For example, for an axisymmetric model consisting 
of a thin stellar disk and flattened dark matter halo, an appropriate choice for F is a thin circular tube in the disk 
plane with (large) radius = 8.5 kpc, the distance between the Sun and the Galactic center. Imagine that V is filled 
with hypothetical observers capable of making local measurements of the dark matter particles. For a given clump, 
we then have / = V~^dV / dp, where dV/V is the the fraction of observers who measure the dark matter density of 
the clump to be between p and p-\- dp. We will estimate / by using N-body simulations to follow the orbital evolution 



and disruption of accreted clumps in the evolving Galactic halo, as described in Section IV 



We adopt a spherically symmetric model for the Galaxy, in which case the remaining parameters in Eq. (2.1) can 
be simplified considerably. Deviations from spherical symmetry alter the orbits of individual clumps but since we are 
interested in properties of an ensemble of clumps (e.g., number of and density within streams) the results assuming a 
spherical halo should provide an adequate approximati on. F or spherical models, there is a one-to-one relation between 
a clump's turnaround radius rta and zta (see Section [II B). The sole remaining parameter required to fully specify 
the orbit of the clump is the specific angular momentum J at turnaround. For a spherically symmetric Galaxy model, 
the local volume V can be replaced by a thin spherical shell of surface area S and radius r^. Eq. (2.1) then takes the 
form 



dP If ,, .dS /r, 

A/ (zta, J, m) — dzta dJ dm . (2.2) 



dp A-Krl J ' ^ dp 



Despite the symmetry, evaluation of Eq. (2.2) appears daunting, since one must sample the space of initial conditions 
(zta, J) for all clump masses m and in each case determine dS/dp from the output of a separate N-body simulation run 
to the present epoch t = to. This difficulty is alleviated by noting that the present state for the large space of orbital 
initial conditions can be sampled by considering a smaller set of orbits at various times. In Figure |l|, for example, 
instead of following the different orbits a,b,c and evaluating them at t = to, it is possible to follow a single orbit 
and evaluate it at three different times such that the dynamical states at a' , b' , c' correspond closely to those above. 
Technically, it can be justified as follows: Let z„ and r„ denote the turnaround redshift and associated turnaround 
radius of a clump that reaches perigee today on its nth passage through the inner parts of the Galaxy. The simulations 
are performed with initial conditions selected from the set (r„, z„), with the provision that dS/dp is evaluated by 
performing an integral over time t taken from the time tn of apogee before the nth passage through the inner Galaxy 
(tn < to) to the time tn+i of apogee after the nth passage (tn+i > to). As illustrated in Figure |l], this essentially 
corresponds to a change in Eq. (2.2) from an integration over zta to an integration over t with a sum over n. 



In principle, Zn depends on the angular momentum J through the usual orbit equations. Howeve r, J must be 
relatively small, since r„ ^ for 1 < n < 4, and therefore the dependence of Zn on J is negligible. Eq. ( |2.2| ) can then 
be written 
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dP A 1 dzt 



d In p Hn dt 

n—l 

where 



/„(J, m) / N {zn , J,m) dJ dm (2-3) 



/„(J,m) = J!^ [ ^dt, (2.4) 



47rr^ y d In p 

and where, for convenience, the probabihty distribution function is now defined with respect to hi p rather than p. 
The quantity /„ is dimensionless and wiU be used extensively in the discussion that follows. To evaluate Eq. ( ^.3| ), 
we select representative clumps characterized by to, J, and n, and follow their evolution via N-body simulations in 
a time-dependent model gravitational potential for the Galaxy. Hypothetical observers located on S measure the 
density as the different clumps pass by, allowing one to determine /„ numerically. 

The response of a dark matter detector to particles in a clump or stream depends on their velocity distribution 
as well as their density. We can estimate the velocity dispersion in a stream using Liouville's theorem, which states 
that the density of particles in phase space is conserved. Consider an infalling, initially virialized clump of mass M 
with an initial characteristic density pi and two-dimensional velocity dispersion Ui = (47r/3)^^^ G^/^M^/^p^^. Here, 
Pi is the virial density, which we take to be 200 times the critical density of the Universe at the formation time of 
the clump. Suppose that the density of the final disrupted clump at the detector \s po- By Liouville's theorem, the 
corresponding velocity dispersion will be 



M V" ( pn ^ - 



^ 1/3 

CD — \ 

Pi 

^ 30kms-M^;^^) (^^) (^^) , (2.5) 

ViQioMoy V 0.03 PS y V 0.03 PS y ' ^ ' 

where ps = CSGeVcm"-^ = 5.4 x 10~^^gcm~'^ is the estimated value for the mean (background) density of dark 
matter particles in the solar neighborhood. The fiducial value of O.OSps has been used, since the streams that are 
likely to have the biggest impact in a detection experiment have a present density po ~ 0.03 ps (see Section^. The 
last factor on the right-hand side of H.5 depends on the formation time of the clump, but only weakly. The essential 
point is that the velocity dispersion of the disrupted clump is significantly less than the bulk velocity of the clump 
particles relative to the Earth, which is typically several hundred kms~^ for recently accreted clumps. The clump or 
stream therefore appears as a "cold" , high- velocity distribution of particles. 

III. GALACTIC HALO MODEL 



In this section, various quantities required to evaluate Eq. (2.3) are derived 



A. Growth of the Galactic Halo 



We assume that the growth of dark matter halos proceeds by hierarchical clustering, i.e., that halos are assembled 
from smaller systems through merger and accretion events. The rate at which the Galactic halo accretes matter enters 
the calculation of dP/dp in two ways: (i) directly, through the factor A/", and (ii) indirectly, since the potential well 
of the Galaxy is becoming deeper with time, causing the orbits of individual clumps to contract. 

The extended Press-Schechter formalism [|o[ ^ provides a phenomenological model for the growth of a dark matter 
halo; it is based on the assumption that nonlinear objects undergo spherical collapse, starting from Gaussian initial 
conditions. In this model, if the present mass of the Galactic halo is Mq, then at redshift z the average number of 
progenitors of this halo that have mass between to and to -|- dm is given approximately by [^ , ^ 



dm V2^ (E(to) - ll{Mo)f'^ 



exp 



{5,{z) - 5M) 
2 (E(to) - S(Mo)) 



dm 



dm 



(3.1) 
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where S(M) is the variance in the present density fluctuation field smoothed in a top-hat window of radius R — 
{3M/4:7rpoY^^ , Po is the present mean density of the Universe, and Sc{z) is the ampfitude that a hnear density 
perturbation, extrapolated to the present epoch, must have in order for the associated object to reach turnaround 
by redshift z. In a matter-dominated, Einstein-de Sitter (f2o — 1) Universe, Sc{z) — Scq (1 -I- z), where Sco = Sc{to) = 

0.15 (Stt^)^^'^ ~ 1.44. Note that this differs from the more common value of Scq — 1-69 since the latter is associated 
with the collapse time, not the time of turnaround. In other cosmological models, the time-dependence of 5c differs 
from the expression above, while the value of Sco is relatively insensitive to cosmological parameters (e.g., ||4l[|). For 
definiteness, we assume a cold dark matter model Universe with Qq + = 1, ^^o = 0.3, and h ~ 0.7, where h is 
the present value of the Hubble parameter in units of 100kms~^Mpc^^, and flo and Q\ are the contributions to 
the total mass density, in units of the critical density pcrit, from non-relativistic matter (dark matter and baryons) 
and the cosmological constant. The primordial power spectrum is chosen to have the scale-invariant form expected 
from inflation, Pi{k) ~ fc, and the present spectrum is normalized to agree with the observed cluster abundance, i.e., 
as = 0.90. (This particular value is obtained from the algorithm described in ji^.) This choice of cosmological model 
fixes the expression for 5c{z) and also determines I](M). 



The average growth rate of a dark matter halo can be determined using eq. (3T). Nusser & Sheth ||4^ (see also 
pl| , |45[ |) have developed a sirnple algorithm for generating possible histories of the most massive progenitor, several 
of which are shown in Figure 0. Though there is significant stochasticity in the various histories we assume a simple 
power-law model for the growth in mass of the Galactic halo 

M{t) = Mo (Pi , (3.2) 



to 

where a = 1.5 (solid line in Figure H) for the cosmology that we have used. We adopt a Milky Way halo mass 
Mq = 2 X IO^^Mq, the value derived in from observations of globular clusters and satellite galaxies. The exact 
details of the halo growth do not significantly affect our calculation of dP/dp, since it is the long term deepening 
of the potential well of the Galaxy that is most important, not the short term fiuctuations in the merger rate. In 
addition, we make the reasonable assumption that the angular momentum distribution for the accreting clumps is 



independent of their mass, so that the J-dependence of Af separates out. This assumption, together with Eqs. (3.1) 



and (3.2), yields the following approximate form for Af 



Af {zte.,J,m) ^ aMo (^] ^xiJ,zte.) , (3.3) 
yiQ / dm 

where x( J, zta) is the normalized clump distribution as a function of angular momentum at turnaround 
(/ dJx{ J, Zta) = 1)- In principle, x may be determined from simulations or from a detailed analysis of tidal 
torques in the hierarchical clustering scenario. For the present discussion, we assume a scaling form for namely 
x( J, Zta) = x(.J I Jcacizta)) where Jcirc(-Zta) IS ttie angular momentum for a clump in a circular orbit at the radius r^a- 
In particular, we assume that the clumps are uniformly distributed in J'^. x then takes the form 

2/32 

X( J, Zta)dJ = e ( J,„ax - J) „ JdJ (3.4) 

where Q{x) is the Heaviside function, /5 = Jdrc/ Jrrmx, and Jmax is a model parameter that characterizes the spread 
in angular momentum of the clumps. We adopt a value of /3 = 2, implying that the maximum angular momentum 
that a clump can have is one-half that required for circular orbits. It is trivial to consider different choices of (3 since 
dP/dp scales with 0^ . At most, with our choice of /3, we overestimate dP/dp by a factor of 4. (This would require 
merging clumps to be on circular orbits which is very unlikely.) In fact, (3 could be significantly higher than 2 if 
inf ailing clumps are on predominantly radial orbits. 

Though halos are never in true virial equilibrium, approximate equilibrium is reached provided no major merger 
events have occurred in the recent past. It is then customary to set the turnaround radius equal to twice the virial 
radius, r2oo, where the latter is defined as the radius within which the mean density of the halo is 200pcrit(t)- For the 
ACDM model, the turnaround radius can then be written 



"'^^^^ " ^800^1^1(0) 



1/3 



520 (no{l + zf + l-no] ' (^) kpc . (3.5) 



-1/3 / < \ "/^ 
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It is well known that the Press-Schechter mass distribution used above does not agree precisely with results from 
N-body simulations (e.g., The mass function obtained from N-body simulations generally has more high-mass 

objects and fewer low-mass objects than the Press-Schechter distribution predicts. Correcting this difference would 
have little effect on our final value for dP/ dp. The more massive infalling clumps would create larger tidal tails, 
thereby increasing the probability of encountering one today. On the other hand, there would also be fewer low-mass 
objects falling into the Milky Way. Taken together, these effects should partially cancel, suggesting that that the 
error introduced by using the Press-Schechter distribution instead of the numerically determined distribution should 
be small. Moreover, the Press-Schechter expression for the progenitor distribution at moderate redshift (III.l) is more 
accurate than the Press-Schechter mass function at late times. 



B. Galactic Potential 



Our simulations follow the evolution of an individual clump in a rigid, time-dependent gravitational potential 
designed to represent the Galaxy. We adopt a three-component model for the Galactic potential out to a truncation 
radius, rtrunc, which is determined by the total mass of the system. Beyond the truncation radius, the potential is 
assumed to be Keplerian: 

sphcr 

-I- $disk if r- < 

^trunc /o c\ 

~ 1 -^ + $00 if ?->?-tm„c ^ ■ ' 

The model halo is described by a logarithmic potential, the spheroid by a Hernquist potential |Q, and the disk by 
the spherical analog of a Miyamoto-Nagai potential . While the shapes of dark halos are not well constrained, the 
results of Ibata, et al on the tidal stream associated with Sagittarius suggest that the Milky way halo potential 
is close to spherical. The assumption of a spherical disk potential is less realistic but should not affect our results 
significantly: while a planar disk would cause the orbits of an otherwise spherical model to precess and leave the orbit 
plane, it will not affect the properties upon which our subsequent calculations are most dependent - the number and 
density of the tidal streams. The components of the Galactic potential in our model are thus: 

*haio = ^t;haioln(r-'+a') , (3.7) 
<i>.phe.. - -^M^ , (3.8) 



*disk = 77^ , (3.9) 



where Uhaio, and a are all tim e-de pendent , with M{t) given by Eq.(3.2). We assume that rtrunc and a scale with 
time at the same rate as rta (Eq. (|3.5| )). The time dependence of Whaio is then set by the relation 



= ^iL^ ,3.10, 



The disk and spheroid potentials are assumed to be time- independent. We use values for the present-day parameters 
that differ slightly from those found in Ref.||5|]: 'yhaio(io) = 200kms"\ a(io) = 16.5kpc, M^phcr = 3.4 x IO^OMq, 
h = 0.7kpc, Mdisk = IO^^Mq, and c = 6.5kpc. For these model parameters, the circular speed at is 220kms~^, 
in agreement with the accepted lAU value, and the rotation curve is relatively flat out to rtrunc- Assuming A/q = 
2 X lO^^M©, the current value of the truncation radius is rtrunc = 216 kpc. 

Selected orbits for which the clump reaches perigee at the present time, io = 13.5 Gyr, on its first, second, third, 
or fourth orbits through the Galaxy are shown in the right panel of Figure |^. To construct these orbits, the (rta, J)- 
parameter space is sampled at random, and those initial conditions for which the orbits satisfy r < at t = io si's 
marked in the left panel of Figure ||. As noted above, the decrease in apogee with time is due to the time-dependent 
nature of the Galactic potential. The orbits are followed from turnaround, which occurs at redshifts zi — 0.26, 
Z2 = 0.56, Z3 = 0.79, and 24 = 1.0 respectively. 
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C. Structure of Clumps at Turnaround 



We next turn our attention to the structure of the clumps at turnaround, before they have been subjected to the 
tidal fields of the Galaxy. We treat the clumps as composedpurely of dark matter, i.e., we ignore the dynamical 
effects of baryons in the clumps. Numerical simulations (e.g., Q, suggest that the density profiles of dark matter 
halos have a 'universal' shape characterized by an inner power-law cusp and an density fall-off at large radii. 
These include the NFW profile and that proposed by Moore et al. [Q. However, for convenience, we model the 
infalling clumps as Hernquist spheres E^, for which the density profile is given by 



3 . (3.11) 

Pcrit r (r + a) 

where pcrit is the critical density for closure, a is the scale radius of the halo, and is the characteristic density in 
units of Pcrit- Since the density falls off asymptotically faster than r~^, this model has the practical advantage that the 
total mass is finite, AI — 271 P cr ito!^ , without having to impose a truncation radius. Furthermore, the corresponding 
particle distribution function, f{E), can be expressed analytically for the Hernquist model |Q, while no closed form 
is available for the NFW or Moore profiles (however, see |5^, [sijl). 

The characteristic density and scale length a are determined using the algorithm outlined by Navarro, Frenk, and 
White iQ. Given the virial mass, M200, and the redshift at which the halo is identified (in our case the turnaround 
redshift), both the virial radius, r2oo, and characteristic density, ^c, can be calculated, independent of the halo model 
assumed. The virial mass and virial radius are defined through the relation 

_ 800 3 , , 

^^200 = — T^rjooPcrit ■ (3-12) 

and r2oo can then be used to c alculate the scale length, a, of the Hernquist profile. For the Hernquist model, the 



virial mass, as derived from ( 3.11 ) is 



M200 = 27recPc„ta'7 ; (3.13) 

(?'200 + ay 



when combined with the definition of M200 (Eq. (B.12)) it yields 

ee=^x(l+x)2, (3.14) 

where x = r2Qo/a. Thus, once and r2oo are known, the scale length a can be readily calculated. 

While we implicitly assumed a 1-to-l correspondence between and a, the results of N-body simulations show 
that there is an intrinsic scatter in the £,c — a relationship [Q. This stochasticity will not significantly alter the results 
since these fluctuations are largely erased when we integrate over the ensemble of infalling clumps. 

The details of the density profile for the clumps (e.g., Hernquist v. NFW) should not significantly alter our results. 
As shown in Figure |^, the central profiles have the same behavior, while the NFW profile falls off less rapidly 
at large r. However, it is the intermediate region, where the density is close to the background dark matter density 
in the solar neighborhood (the shaded region in the Figure) that is of greatest interest. The central cores of the 
accreting clumps remain compact while the outer regions are quickly stripped by the tidal field of the Galaxy. Since 
the density fall-off is more gradual for the NFW model the resulting tidal streams would be longer, slightly increasing 
the probability that the solar system is in a low density stream today. Thus, the choice of a Hernquist profile should 
underestimate dP/dp at low clump densities. 



As described above, the calculation of dP/dp requires that we evaluate the quantity /„ in Eq. (2.4) for various 
clump parameters. An approximate form for /„ is found by assuming that the clump is unaffected by the tidal field 
of the Galaxy. In this limit, 

fn = ^° (3 15) 

where Ur,n is the radial velocity of the clump on its n'th pass through the solar neighborhood (see Figure ^). The 
quantity dV/dp is readily calculated from the clump density profi le. The resulting expression for /„ can be compared 
directly with the measurement of /„ in the simulations (Eq. ( |2.4|) ), as we will see in the next section. 
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IV. NUMERICAL SIMULATIONS OF CLUMP EVOLUTION 

We use numerical simulations to study the effects of tidal fields on clumps as they pass through the inner regions 
of the Galaxy. Clumps of various masses, initially described by the Hernquist spheres, are set on orbits such that a 
point particle with the same initial conditions would reach the solar radius today. The clumps are followed as they 



move through the rigid, time-dependent Galaxy potential described in Section [II B 



A. Code Implementation and Density Measurements 

Our simulations use the Barnes-Hut treecode ||55|] modified to incorporate a background potential and to allow the 
calculation of densities in situ at specified time intervals. The density at any point is approximated by taking the 
weighted average over N nearest neighbors: 

N 

p(r0^5]m,W^(r„r,) , (4.1) 
where nij is the mass of the j'th nearest particle, W (r^, rj) is the symmetric smoothing kernel, 

W{r^,r,) = ^^{w{\r,-r,\/H,) + w{\r,-r,\/H,)) , (4.2) 
where Vi = ATrHf/3, Hi is half the distance to the A^'th nearest neighbor to ab observer at r.i, and 

if X < 1 

w{x) ^ { \(2-xf if 1 < X < 2 (4.3) 

otherwise 

This method is used extensively in smooth particle hydrodynamics simulations. 

Clumps are modeled with 10^ particles. The time step At is dynamically adjusted pfl] using 




At^aJ^, (4.4) 

where e is the softening length, Cmax is the maximum acceleration that any particle has for that timestep, and a ~ 0.5 
is a fixed parameter adjusted to optimize the performance of the treecode. The softening length is taken to be l/40th 
of the scale radius of the infalling halo. We have checked that the results do not change significantly when the number 
of particles is increased or the timestep reduced. 

The quantity /„ is calculated by placing fictitious observers on a spherical shell of radius — 8.5 kpc centered 
on the model Galaxy. At regular intervals throughout the simulation, these "observers" record their local density 
using the method described above. 10"* observers are used in order to ensure that no significant stream or clump slips 
through the surface S undetected. 

We note that the effects of dynamical friction on the evolution of the clumps are not included in the code. Dynamical 
friction due to the motion of a clump through the Galactic halo would cause a steady deceleration of the clump in 
the direction of its motion, leading it to spiral into the center of the Galaxy. These effects should be negligible for the 
problem at hand, since the timescale for dynamical friction is generally much longer than the age of the Universe. For 
example, for a clump of mass M ~ IO^Mq travelling at speed Wciump = 500 km s^^ through the halo, with a perigee 
of r = 8.5 kpc, the timescale for the clump to spiral into the Galactic center is at least 



(Cf. Binney and Tremaine, Section 7.1). This is a conservative underestimate of the friction timescale because it 
assumes the clump spends all its time at perigee and ignores the spatial extent of the clump. A clump more massive 
than this which penetrates the solar orbit could in principle suffer significant dynamical friction. However, even in 
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this case, we expect the chimp to be ripped apart by the halo tidal field before it is appreciably slowed by dynamical 
friction; as a result, the effective clump mass (as far as friction is concerned) is reduced, again rendering friction 
unimportant. To see this, a crude estimate of the ratio of the tidal and dynamical friction forces gives 



lump 



where 7 is the ratio of the clump density to the local smooth halo density. Furthermore, if we were to include 
dynamical friction, it would cause more clumps to pass through the solar neighborhood (on their way to the Galactic 
center), increasing the probability of detecting a clump. 



V. RESULTS 

Figure || provides snapshots of lO^M© clumps that reach the solar neighborhood today on their first, second, third 
and fourth orbits through the Galaxy. The effects of tidal fields on the clumps are clearly evident. During the first 
couple of orbits, a clump is stretched into a long tidal stream which, by the third and fourth orbits, wraps around the 
Galaxy. However, the high density clump cores survive for many dynamical times. 

Figure ^ illustrates what our sphere of hypothetical observers see as a clump passes through the inner part of the 
Galaxy for the first time. There are two distinct regions of high measured density where the clump enters and exits 
the r = Vs sphere. Thus, at this instant, a small fraction of observers measure a high density of high-velocity clump 
particles while most of the observers measure a relatively low density of clump particles. 

Figure ^ shows the measured values of /„ as a function of p for IO^Mq clumps that reach the inner region of the 
Galaxy by the present epoch on their first, second, third, or fourth orbit. The analytic result (Eq. ( 3.15| )) is shown 



for comparison. For clumps on their first passage through the inner region of the Galaxy, the numerical and analytic 
results for /„ as a function of p have roughly the same form. This is to be expected since tidal fields have only 
just started to disrupt the clump by this stage in its evolution. The analytic result underestimates the measured 
values by a factor ~ 4. This is easily understood as follows: In evaluating Eq. ( 3.15D we set Vr = 550 kms~^, the 



characteristic velocity of the clump at the solar radius. In fact, we should use the component of the clump velocity 
normal to the sphere of observers which, for orbits whose turnaround radius is comparable to r^, is considerably less 
than 550 kms~^. 

Tidal effects distort the shape of the /„ vs. p curve. The general trend is to increase the probability at low densities 
and decrease the probability at high densities, the transition occurring at p = 10~^ pB- Note that for very high 
densities, /„ approaches the value expected from the analytic results. These regions are at the centers of the clumps 
where the densities are so high that they are not tidally disrupted. 

Our result for P{pd > p), the probability that the Earth is passing through a stream of density greater than p is 
shown in Figure ^. This figure indicates that there is a high probability for the Earth to be passing through a clump 
or stream with density 3% of the mean local halo density pB- The probability drops precipitously with density to 
values below 0.01 for streams that would contribute a density comparable to that of the background. 

The separate contributions to dP/dp from objects on their first through fourth passages through the Galaxy are 
shown in Figure ^. The largest contributions come from objects that have made several passes through the Galaxy. 
This is primarily due to the increase in accretion rate with turnaround redshift. Additionally, clumps are stripped 
of particles with each passage through the Galactic center so that the tidal stream of, say a "fourth-pass clump" 
wraps several times around the Galaxy. Thus, the contributions to the density can come from what are essentially 
independent phase space streams. One might be concerned that our assumption of a spherically symmetric Galaxy 
potential artificially forces these streams to lie in a single plane, enhancing dP/dp at large p. However, we are confident 
that our cut-off at fourth-pass clumps is sufficiently conservative that this should not be a big effect. 

Finally, we turn to the velocity distribution of particles as a clump or stream passes through the observation shell 
S. The velocity space coordinates of particles within 0.5 kpc of S for the four snapshots in Figure I (10^ Mq clumps 
that have made 1 — 4 passes through the Galaxy by the present epoch) are shown in Figure |l^. On its first pass 
through the inner part of the Galaxy, the clump is tidally disrupted with particles spreading out along the orbit of 
the clump center-of-mass. In velocity space, particles have a fairly tight distribution in speed but are spread out in 
direction as we expect from Liouville's theorem. By the third pass through the Galactic center, several streams are 
apparent and the velocity distribution becomes rather complicated. 

More relevant for our purposes is the velocity space distribution as determined by a given observer. For illustrative 
purposes, we choose observer 'O' located at the instantaneous position of the center of mass for the clump in Figure ^ 
that has made three passes through S. Figure |ll| shows the velocity space distribution for particles within 0.5 kpc of 
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O. The distribution is anisotropic with an rms dispersion along the direction of motion of 19knis^^ and dispersion 
along the two transverse directions of 34kms~^ and 53 kms~^, again consistent with what one expects from Liouville's 
theorem. A more detailed discussion of the evolution of a satellite's velocity space distribution and the connection 
with Liouville's theorem can be found in pal. 



VI. EXPERIMENTAL SIGNATURES OF ACCRETED CLUMPS AND STREAMS 



A. WIMP Search Experiments 



At present, over twenty groups around the world have deployed or are in the process of building terrestrial detectors 
designed to search for WIMPs ||5^. These experiments attempt to measure the energy deposited when a WIMP 
interacts with a nucleus in the detector. A WIMP of mass that scatters elastically with a nucleus of mass tojv 
deposits an energy Q = {mlv^ /rriM) (1 — cosO*), where = mf^m^j {iriN + m^) is the reduced mass, v is the speed 
of the WIMP relative to the nucleon, and 9* is the scattering angle in the center-of-mass frame. The differential 
detection rate (per unit detector mass) can be written (e.g., Q): 



where tJo is the scattering cross-section, pL is the local WIMP density, G{Q) is a form factor for the WIMP-nucleon 

1 /2 

interaction, and fmin = [QniN /'^rn'^) ■ Here F{v) is the normalized distribution of WIMP speeds (/ F(v)dv = 1) 
in the rest frame of the detector, obtained by integrating the three-dimensional velocity distribution /(v) over angles. 

In keeping with the discussion of the previous section, we assume that the local distribution of WIMPs can be split 
into a smooth background and a single coherent clump or stream, so that 



dR _ dRs dRc , , 

dQ ~ dQ dQ ^ ' ' 

- '""^^ G^Q){Tb{Q)+Tc{Q)) , (6.3) 



where the subscripts B and C refer to the background and clump respectively, and ps denotes the smooth local halo 
density in the absence of the clump. The remaining terms are given by 



Tb{Q) 



Fb{v 



dv 



(6.4) 



and 



TciQ) = ^ r Fciv)^-^ . (6.5) 

A standard model for the halo background velocity distribution is the Maxwellian distribution, assumed to be 
isotropic in the rest frame of the Galaxy: 



/b(v') 



^ / 2 



(6.6) 



where v' is the WIMP velocity in the rest frame of the Galaxy, and f s is the two-dimensional halo velocity dispersion, 
corresponding approximately to the circular velocity at the position of Sun in the Galaxy, vb — 220 km s~^. Let 
V^;^ be the velocity of the Earth relative to the Galaxy frame (assumed to be non-rotating), so that v' = v + V^s- 
Transforming the distribution (VI. 6) into the Earth frame (e.g., |5^) and substituting into (VI. 4) yields 
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In what follows, we use a slightly more complicated expression that assumes a velocity-space cut-off at the escape 
speed of the Galaxy with the value Vcsc = 575 kms""'^ (e-g-, |5^)- If we assume that the velocity distribution of the 
WIMPs in a clump is Maxwellian and isotropic in the rest frame of the clump, with internal velocity dispersion vc, 
then 

TciQ) = (eJYE£±IlIlill.) +crff ^"^-"- )) (6.8) 

^PbVec \ \ Vc J \ J J 

where Vec is the relative velocity of the Earth and the clump. 

It is straightforward to understand heuristically the effect of a WIMP clump on the nuclear recoil energy spectrum 
in a dark matter detector. The key observation is that the clump internal velocity dispersion vc ~ 50 kms~^, is small 
compared to Vec ~ 525 km s^^. In the limit vc/Vec ~^ 0, the clump particle distribution function in the frame 
of the Earth is given by Fc{v) cx S{v — Vec)', from (VI. 1), this implies that the detection rate is constant, with an 
amplitude proportional to I/Vec, if Vec > Wmin(Q) and zero otherwise. That is, the detection rate is constant at 
small Q, and zero beyond some maximum value Qmax = ^m^V^c/^N- The value of Qmax varies seasonally, as the 
Earth's velocity relative to the clump changes. Including the non-zero clump velocity dispersion vc turns this sharp 
cutoff in the detection rate into a smooth shoulder, with a width proportional to vc- 

As discussed above, the distribution of particles in velocity space for a tidally disrupted clump is anisotropic. In 
particular, the dispersion is a factor of 2 — 3 times smaller along the orbit of the clump (defined by ~Vcb, the velocity 
of the clump in the Galaxy frame) as compared with the transverse directions (cf. Figure 11). Fc(v) is calculated in 
the Earth frame of reference (in Figure 10, the distribution of v about a point displaced from the origin by Veb)- The 
width of the peak of Fc{v) therefore depends on the angle between Veb and Vcs- When these velocities are colinear, 
the distribution is relatively narrow. When they are close to perpendicular, the distribution is somewhat broader. In 
either case, the dispersion in v is small compared with the dispersion of particles in the background. Thus, velocity 
space anisotropy alters the detailed shape of the feature in a recoil spectrum but not its general characteristics. While 
it is straighforward to include velocity space anisotropy in the calculations that follow, doing so increases the number 
of parameters and complicates the discussion. Therefore we work under the assumption that the clump velocity 
dispersion is isotropic. 

The recoil energy spectra for Vec = 401, 569, and 698 kms""'^ are shown in Figure 12. These velocities correspond 
to a clump moving with a speed of Vcb — 525 kms^^ relative to the Galaxy at an angle with respect to Veb of 
45°, 90°, 135°, respectively. For Veb we use the Sun's velocity relative to the Galaxy (seasonal modulation will 
be discussed below). In this example, we have assumed a clump with local density equal to 3% of the local halo 
density, i.e., pc = 9 x 10"'^ GeVcni^-^ at the position of the Earth, which Fig. 8 shows to be a common occurence. 
dR/dQ is given in units of Events per kg of detector per keV per day with assumed parameters ctq = 4 x 10~^^cm^, 
m-^ = 50 GeV and niN = 68.5 GeV corresponding to a ^'^Ge detector. In addition, we have set the nuclear form factor 
G{Q) = 1. As anticipated, the contribution from the clump is well-approximated by a step function in Q with a 
more pronounced contribution from the clump at larger values of Vec- Though the number of events per energy bin 
decreases as 1/Vec, the total number of events increases as Vec- A positive feature of the clump signal, in terms of 
detect ability, is that it dominates over the background halo signal at high recoil energy, where radioactive backgrounds 
are generally small. A negative feature is that the absolute detection rate associated with a clump is not likely to be 
large (which is why it is subdominant at lower energies), so a large total detector mass would be necessary. 

Figure 13 shows recoil spectra for various WIMP masses. The background spectra exhibit the usual dependence 
with m^: spectra for smaller fall off faster with Q. Likewise, the clump feature moves to higher energies as 
is increased. 

While the hierarchical model gives statistical predictions for the clump density, velocity, and velocity dispersion, the 
seasonal variation of the signal depends on the relative velocity of the clump and the Sun, which we cannot predict a 
priori. Figure 14 shows several choices for the clump or stream velocity vector. The direction of motion of the Sun 
around the galaxy and the orbit of the Earth around the Sun (the plane of which is perpendicular to the page) are 
also shown. The Earth is at the upper left point of its orbit on about June 2 and at the lower right six months later. 
For directions A or B, the annual modulation of the clump signal, that is, the amplitude of the yearly variation in the 
maximum recoil energy Qmax, is maximal, of order 25% peak to peak; for streams C and D, perpendicular to Earth's 
orbit, there is no seasonal modulation of the clump signal. For stream A, Qmax reaches a minimum on June 2 and a 
maximum six months later; for stream B, the situation is reversed. Note that the four directions shown here all lie in 
the plane defined by the Sun's velocity and the orbital axis of the Earth around the Sun. In such cases, the annual 
modulation of the clump signal is in phase (or anti-phase) with the annual modulation of the background halo (which 
is set by the relative velocity of the Earth with respect to the Galactic halo), with a peak or trough on June 2. We 
emphasize that this is a special case: in general, the clump orbit will not lie in this plane and will therefore have a 
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different phase from that of the smooth background halo. Consequently, the phase of the clump signal modulation 
constrains the clump direction of motion. 

Figure 15 shows the seasonal modulation for example B of Figure 14. The two solid lines show the shifts in the 
spectrum due to seasonal variation of the clump signal, while the seasonal modulation of the background signal [s^ 
is shown by the two sets of dotted curves. The fractional modulation signal. 



idR/dQ),^^, + {dR/dQ)^^^ 

is shown in Figure 16 for both B (clump feature and background in phase) and A (feature and background out of 
phase) . 

If the recoil spectral feature due to the clump can be detected, then it should be possible in principle to determine 
the physical parameters of the dark matter stream: its density, direction of origin, velocity, and velocity dispersion. 
The maximum recoil energy associated with the clump Qmax determines the product rrirVEc- The background halo 
signal itself determines rrirVB', if vb is sufficiently reliably determined from the galaxy rotation curve, then Vec can 
be extracted. The clump velocity dispersion vc is then determined from the width of the clump recoil shoulder near 
Q = Qmaxi and the clump density pc from the relative amplitude of the clump and background halo signals. Finally, 
the phase of the clump modulation signal constrains the plane of the clump orbit, while the amplitude of the clump 
modulation signal, along with the parameters above, constrains the direction of the clump velocity vector in that 
plane. In practice, of course, given finite detector energy resolution, the presence of radioactive backgrounds, low 
counting rates, and possible confusion due to more than one clump, the dynamical information that will be extracted 
is likely to be more limited. 

A new generation of WIMP detectors with directional sensitivity (e.g., DRIFT [^) is now under development. 
The ability to measure the nuclear recoil direction would provide a new avenue for constraining the clump speed 
and direction relative to the Sun. Fig. |l^ shows the angular recoil spectra for two special cases: clump travelling 
through the Halo in the same direction as the Sun and clump travelling in the opposite direction from the Sun. In 
each case, we assume that the clump speed in the Galactic rest frame is 525 km s^^. We choose a clump dispersion 
Vc — 25kms~^ which is a factor of two smaller than what was used in the previous discussion. (Recall that the 
velocity dispersion is smaller along the direction of motion of the clump than in the two transverse directions. Since 
we are considering the special case where clump and solar system velocities are colinear, we select a dispersion at the 
low end of the expected range.) When the WIMP stream is travelling in the same direction as the Sun there is an 
increased probability of forward scattering. When the stream travels against the sun's motion, rearward scattering is 
increased. Note that there is an appreciable signal when clump and Sun are moving in the same direction whereas 
the feature in the energy spectrum for this case is barely detectable. 



B. Axion Search Experiments 



Axions are pseudo-scalar particles that arise in theories designed to solve the Strong CP problem (e.g., |]60| ). 
Although in most models these particles are extremely light, in the early Universe, they are non-relativistic and 
therefore behave like other cold dark matter candidates. The direct detection scheme for these particles relies on their 
coupling to photons (077) ||6^, In a magnetic field, the axion can convert into a microwave photon with an energy 
equal to the axion's total (rest mass plus kinetic) energy; this signal can be detected with a suitably tuned microwave 
cavity. In this respect, axions would provide a more direct picture of the distribution function of dark matter. The 
axion energy spectrum in Figure |l^ shows that significant features can be introduced into the spectrum (thick lines) 
which are not present in the pure Maxwellian background model (thin lines). Here, we have assumed the same 
parameters as in our discussion of WIMP detection via nuclear recoil {vc ~ 50kms^^, Vec — 569 km s^^, pc/pB = 
0.03). The feature appears as a platueau at high velocities, right where the background is dropping rapidly. Moreover, 
if the clump dispersion is smaller, the feature appears as a distinct secondary peak. This is illustrated by the dotted 
curve where we assume vc — 25kms~^. 

By following the motion of the secondary peak throughout the course of a year, considerable additional information 
about the dynamics of the stream can be determined. The range of velocities over which the peak moves corresponds 
to the angle of the stream with respect to the plane of the Earth's orbit. For example, a stream perpendicular to the 
Earth's orbit would show no variation over the course of a year while one parallel would show the maximum variation. 
The phase of the seasonal modulation in the position of the peak can be used to determine the azimuthal direction 
of the stream as well. This would give a complete dynamical description of the stream. Since the axion detection 
spectrum gives direct access to the axions' energy distribution, it is likely easier to determine the stream parameters 
than for the case of WIMPs. 
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VII. CONCLUSION 

Dark matter detection experiments may provide a new window to the Universe. If successful, the results would 
have dramatic implications for the field of particle physics, leading to physics beyond the Standard Model. The 
implications for cosmology would be equally important since dark matter represents ^ 30% of the energy density of 
the Universe and ~ 90% of the mass in gravitationally bound structures. The results presented here suggest that 
dark matter search experiments can also provide clues to the formation history of the Galactic halo. 

Most analyses of direct detection experiments take, as a starting point, the assumption that the velocity-space 
distribution of dark matter particles in the solar neighborhood is smooth and approximately Maxwellian. However, 
N-body simulations indicate that CDM halos have a high degree of substructure due to undigested lumps and tidal 
streams. To what extent can this substructure affect the results of dark matter detection experiments? Current N- 
body simulations do not have the resolution to answer this question . It is for this reason that we have developed 
a hybrid approach based on the extended Press-Schecter formalism and N-body simulations of individual clumps 
evolving in a background potential. 

Our main results are summarized in Figure ^. Taken at face value, they imply that there is a high probability 
for the local density of dark matter particles to be enhanced by a clump or stream with a density equal to ~ 3% 
of the background. The probability drops precipitously as a function of clump density and is negligible for densities 
comparable to that of the background. 

A dynamically 'cold' dark matter stream with a density ^ 3% of the background can produce an observable feature 
in the energy or angular spectra of a terrestrial search experiment. For recently accreted clumps and their associated 
tidal streams, the anticipated features appear at high energies precisely where the signal from the background drops 
dramatically. In addition, the seasonal modulation of the stream's signal will be different from that of the background. 
In principle, energy and angular spectra would allow one to extract various clump parameters such as its internal 
velocity dispersions, and velocity (speed and direction) relative to the Earth. We conclude by noting that this model 
of halo substructure could yield other observable detection signatures that would be interesting to explore — for 
example, in events arising from WIMP annihilations in the cores of accreted clumps. 
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FIG. 1: Schematic illustration of the change in variables from zta to t used in evaluating dP/dp. Top: six orbits that reach 
different radii at t = to. Bottom: two orbits used to approximate the contributions bracketed by orbits a — c and d — / of the 
top panel. 
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FIG. 2: Three realizations of tlie growth of a 2 x 10 Mq dark halo. The heavy straight line is the analytic fitting formula 
M{t) = Mo{t/tof-^ used to calculate dP/dp. 
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FIG. 3: Initial condition parameter space and corresponding orbits. Left: rta — J parameter space for orbits that are near the 
surface S aX, t = to- The strips of points at rta ~ 380, 285, 230, and 195 kpc correspond to particles that reach S after 1, 2, 3 
or 4 orbits through the Galaxy respectively. The large dots are the initial conditions used in the numerical simulations. Right: 
Distance from the Galactic center as a function of time for the orbits shown in the left panel. 
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FIG. 4: Density profiles for Hernquist (solid curve) and NFW (dashed curve) models assuming a virial mass, IO^Mq. The 
shaded region corresponds to densities at the detector of p = ^ ps — IO^^/Ob where ps = 6 x lO^'pcrit is the density of dark 
matter particles in the background. 



FIG. 5: Snapshots of IO^Mq clumps whose center of mass is near the surface 5 at t = to on their first (upper left) through 
fourth (lower right) orbits through the Galaxy. The surface S is represented by the sphere in the center of each panel. Darker 
regions correspond to higher densities. The symbols used for the particles are elongated in the direction of their velocity. The 
'O' in the lower left panel signifies the position of an observer used later in the discussion. 
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FIG. 6: Sphere of observers and the density they measure as a stream passes through. Darker regions of the sphere correspond 
to higher measured densities. 




FIG. 7: The quantity /„ as a function of density for 10^ Mq clumps that reach S at t = to on their first through fourth orbits 
through the Galaxy. Solid lines show the analytic results assuming no tidal disruption while the symbols show the results 
obtained from the numerical simulations. The two types of symbols correspond to two different choices of orbital angular 
momentum shown in Figure — crosses are the higher angular momentum, triangles are the lower angular momentum. 
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FIG. 8: Probability for the Earth to be passing through a stream with local density greater than p. p is given in units of the 
background density, ps = 0.3 GeVcm"^ = 8 x 10"^ Mopc"^. 




FIG. 9: Contribution to dP/dlnp from objects on their first through fourth pass through 5* (lower line through to upper line 
respectively) . 



21 



Pass 2 



i 



-500 Q 500 -500 500 




-500 500 -500 500 

V, [km/a] V, [kn-i/s] 



FIG. 10; Distribution in velocity space of particles within 0.5 kpc of the observer shell S for clumps on their first through fourth 
orbit through the Galaxy. Particle distributions are the same as in Figure ^ Velocities are projected onto the orbital (xy) 
plane. 
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FIG. 11: Three projections of the velocity space distribution of particles within 0.5 kpc of the observer 'O' that is indicated in 
the lower left panel of Figure pi Velocities are given in kms^^. 
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FIG. 12: Theoretical differential event rate per kg per keV per day as a function of recoil energy for different clump speeds 
as measured in the rest frame of the Earth. The clump speeds are, from left to right, 401, 569, and 698kms~^ corresponding 
to a clump moving through the Galaxy with Vca = 525 km s"'^ at angles relative to "Veb of 45°, 90° and 135° respectively. 
The dotted, dashed, and solid lines are the detection rates due to the background, the clump and the clump plus background 
respectively. Parameters chosen for the calculation are pc/pB = 0.03, ao — A x 10^^''cm^, G{Q) = 1, — 50 GeV, and 
mjv = 68.5 GeV. 
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FIG. 13: Theoretical differential event rate for WIMP masses m-^ — 25 GeV (dashed line), — 50 GeV (solid line), and 
rrix = 75 GeV (dotted line) . Other parameters are the same as in Figure |l3. 
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FIG. 14: Various stream directions as viewed in the rest frame of the solar system. Also shown is the Sun and its direction 
of motion around the Galactic center and the Earth and its orbit around the Sun. Spectral features due to Streams A and B 
will have the maximum seasonal modulation while Streams C and D will show no seasonal modulation. The modulation of the 
feature due to Stream B (A) will be in phase (180° out of phase) with that of the background. 
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FIG. 15: Theoretical event rate vs. deposited energy for June (red lines) and December (blue lines). Dotted lines are the 
signal from the background; solid lines are the si gna l from background plus 3% clump. Top panel: modulation of stream and 
background are in phase (stream B from Figure Bottom panel: modulation of stream and background are 180° out of 
phase (stream A from Figure h4) . 
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FIG. 16: Fractional modulation signal J-{Q) as defined in Eq. (|6^) vs. Q The solid line shows the result for the background. 
The dashed line shows the result for background plus stream B of Figure H (modulation signal of stream in phase with that of 
background). The dotted line shows the result for background plus stream A (modulation signal of stream 180° out of phase 
with that of background). 
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FIG. 17: Relative event rate as a function of recoil angle for a Maxwellian background (dotted line), a 3% clump (dashed line), 
and clump plus background (solid line). 9 = Q corresponds to a recoil event in the direction opposite to that of the Sun's 
motion around the Galaxy. Velocity of clump and Sun are assumed to be colinear. Left: Clump moving in the same direction 
as the Sun. Right: Clump moving opposite to the Sun. Clump speed in the rest frame of the Galaxy is 525kms~^. Other 
parameters are the same as in the previous Figures. 



26 




FIG. 18: Speed distribution for Maxwellian background (light lines), and clump plus background (heavy lines). Shown are 
distribution functions for June (solid lines) and December (dashed lines) assuming seasonal modulation of clump and background 
are in phase (Stream B in Figure and a clump velocity dispersion vc = 50kms~^. Also shown is the distribution function 
for June assuming a velocity dispersion of 25kms~^ (dotted). A clump speed relative to the Sun of 569kms~^ was chosen. 



